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Current data exclude cosmic strings as the primary source of primordial density fluctuations. 
However, in a wide class of inflationary models, strings can form at later stages of inflation and 
have potentially detectable observational signatures. We study the constraints from WMAP and 
SDSS data on the fraction of primordial fluctuations sourced by local cosmic strings. The Bayesian 
analysis presented in this brief report is restricted to the minimal number of parameters. Yet it is 
useful for two reasons. It confirms the results of 0] using an alternative statistical method. Secondly, 
it justifies the more costly multi-parameter analysis. Already, varying only three parameters - the 
spectral index and the amplitudes of the adiabatic and string contributions - we find that the upper 
bound on the cosmic string contribution is of order 10%. We expect that the full multi-parameter 
study, currently underway Q, will likely loosen this bound. 

PACS numbers: 98.80.Cq 



The idea that inflationary cosmology might lead to cos- 
mic string formation is not new [3j and has received new 
impetus from brane inflation scenarios suggested by su- 
perstring theory 0, H H S S H E3 . In these models, 
inflation can arise during the collisions of branes that 
coalesce to form, ultimately, the brane on which we live 
P. HlLll3 ,IT^. It is possible to build brane inflation mod- 
els that predict adiabatic temperature and dark matter 
fluctuations capable of reproducing all currently available 
observations. A seemingly unavoidable outcome of brane 
inflation, however, is the production of a network of lo- 
cal cosmic strings 0, @, whose effects on cosmological 
observables range from negligible to substantial, depend- 
ing on the specific brane inflationary scenario @, Q • As 
the precision of cosmological observations increases, we 
might hope to be able to distinguish among numerous 
presently viable models of inflation by the properties of 
the cosmic strings they predict. 

The fluctuations resulting from brane inflation are ex- 
pected to be an incoherent superposition of contributions 
from adiabatic perturbations initiated by curvature fluc- 
tuations and active perturbations induced by the decay- 
ing cosmic string network. The resulting CMB tempera- 
ture power spectrum can be written as 

Cl = M/Cfiiabatic + acjrtrfni* 

where W and B are weighting factors. Analogous ex- 
pressions hold for matter density and polarization power 
spectra. In Eq. QJ, the weight factors W and B deter- 
mine the relative importance of the adiabatic and cosmic 
string contributions. We choose the weight factor W so 
that W = 1 for the WMAP's best fit power law model 
[TH with no cosmic strings. 

The scaling properties of the cosmic string networks 
in brane inflationary scenarios are different from those of 



the more familiar 3+1 D field theoretical strings, since 
intercommutation probabilities are smaller as a conse- 
quence of the existence of extra dimensions Hi. As 
argued in observational constraints on the ampli- 
tude of string-induced CMB perturbations place limits on 
Gfi/y/X, where /i is the string tension and A < 1 is a di- 
mensionless measure of the intercommutation rate. The 
constraint reported in based on WMAP M, EH 
and 2dF data H3, was G[i < 1.3 x 10- 6 y/BX/0A, where 
B, defined by Eq. JIJ, measures the importance of pertur- 
bations induced by cosmic strings. It was argued that, 
conservatively, the currently available data still permit 
B < 0.1. This bound on the parameter B was found 
by minimizing the x 2 with respect to variations in the 
parameters B, W, and n s , the spectral index of the adi- 
abatic fluctuations power spectrum. The uncertainty in 
the minimum \ 2 was translated into an upper bound on 
B. This statistical procedure was not rigorous, although 
the derived bound should be reasonably secure. 

In this paper we re-calculate the bound on B, vary- 
ing the same parameters as in Q, but using a rigorous 
Bayesian likelihood analysis. Except for the statistical 
methods used, and for fitting to the SDSS instead of the 
2dFGRS power spectrum, our setup is identical to the 
one used in Q and the interested reader can find most 
of the details there. Here we only mention that C ; s rings 
were computed using a modified version of CMB- 
FAST 19] which uses as its active source the energy- 
momentum tensor components calculated using the lo- 
cal string network model originally suggested in |20j and 
further improved in [l8L l2lT ] . C adlabatlc were computed 
using a recent of version of CMB FAST. As in [l|, we 
vary the value of the spectral index n s , while holding the 
following parameters fixed at their best-fit values as de- 
termined by WMAP [3: n m h 2 = 0.14, Q b h 2 = 0.024, 
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f?A = 1 — — ^b, h = 0.72, r = 0.166. Our rationale 
is that cosmic strings are only a small correction to the 
WMAP best-fit. Elsewhere y], we shall present a more 
complete study in which we also allow fl m h 2 , 0;,/i 2 , t and 
h to vary as well. 

According to Bayes's theorem, the posterior distribu- 
tion for W , B and n s is given by 



P(W, -B, n s \M)C(D\W, B, n s ,M) 
P(D\M) 



P(W,B,n s \D,M) 

r\u\m ) 

(2) 

where C{D\ W, B, n s , M) is the likelihood, M denotes 
our model and D stands for data. We use WMAP's 
CMB temperature anisotropy power spectrum (TT) and 
the temperature-polarization cross-correlation (TE) to- 
gether with the power spectrum (P(k)) from the SDSS 
experiment. For TT and TE we have used the likelihood 
function provided by the WMAP team [lg, [22J and for 
P(k) th e SDSS likelihood code provided by Max Tegmark 
[HlUl. The combined TT+TE + P(k) likelihood is sim- 
ply the product of individual likelihoods. 

We choose fiat priors, P(W,B,n s ) for parameters W, 
B and n s with ranges W € [0.95, 1.05], B € [0.0,0.105], 
and n s £ [0.955,1.02]. We then survey this parameter 
space using a 73 x 65 x 27 grid for W, B and n s , respec- 
tively. At each grid point we calculate four predicted 
spectra from our model: Cf T , Cf E , C EE , and P(k). 
The SDSS power spectrum is determined up to an over- 
all bias factor, so the bias has to be included into our 
analysis as a fourth parameter. The bias and the effect 
of the spectral distortions on P(k) can be parametrized 
via an overall multiplicative factor, which we vary over 
the range [1.50,1.86] - a range large enough to include 
the maximum likelihood bias value for every grid point - 
and over which we marginalize. After compiling the full 
likelihood grid, we marginalize over W and n s separately 
and plot the resulting likelihood contours. In Figs.Q] an d 
Elwe show the contours in the (W, B) and (n s , B) planes 
respectively. Finally, we marginalize over both W and 
n s to find the posterior for B, shown in Fig. [21 In a sim- 
ilar manner, we have also generated the posteriors for 
W (Fig. 0} and n s (Fig. EJ. The posterior results for 
W can be interpreted in light of Fig. 0s results for B, 
especially with the help of Fig. ^ We see that there is 
something of a degeneracy between W and B, so that 
an increase in B is compensated by a decrease in W. 
Nevertheless, our determination, W = 0.985 ±0.025, has 
a much smaller range than that found by WMAP with 
B = 0: W = 1.0 ± 0.11 [H . For the spectral index, 
we fin d ru = 0.995 ± 0.0165. For comparison, Spergel 
et al. [15| find n s = 0.99 ± 0.04 but vary many more 
parameters. We note, from Fig. |2J that the approximate 
Aln(£) = 1 range for n s is » 0.995 ± 0.005 at B = 0. 
We therefore suspect that the relatively narrow range we 
find for n s could broaden when we vary Q m , fib, t and h 
in addition to W, B, and n s . 

Our principal result is Fig. [31 the posterior for our 
cosmic string weighting parameter, B. We note that, de- 
spite our use of the unadjusted WMAP best-fit param- 
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FIG. 1: A(ln(£)) = 1, 2, 3, and 4 contours after the likelihood 
grid has been marginalized over n 3 . The maximum likelihood 
model is marked. 




Spectral Index n a 

FIG. 2: A(ln(£)) = 1, 2, 3, and 4 contours after the likelihood 
grid has been marginalized over W. The maximum likelihood 
model is marked. 
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FIG. 3: This plot shows P(B\D, M) after marginalization 
over both W and n a . The three dots indicate the regions con- 
taining 90, 95, and 99 per cent of the probability, respectively 
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FIG. 4: This plot shows P(W\D,M) after marginalization 
over both B and n 3 . The dots indicate the regions containing 
90, 95, and 99 per cent of the probability, respectively 



eters for all cosmological factors save n s , our additional 
cosmic string degree of freedom still improves the fit to 
the data, with the likelihood maximum likelihood falling 
at B « 0.04. However, at the 90% confidence level, our 
results are also roughly consistent with B = 0, while the 
upper bound at 99% confidence is B < 0.09, quite close 
to what was found much less rigorously in This up- 
per bound might still permit observable signatures such 
as pulsar timing |2E| . gravitational wave bursts [2(| , and 
distinctive double lensing events [13, HU . That the peak 
for B is not at zero, however, underscores the necessity of 
varying the other cosmological parameters, r, h, VL mi and 
Of, as well. We have begun work on this multi-parameter 
analysis using the Markov Chain Monte Carlo method 
0. 
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FIG. 5: This plot shows P(n s \D, M) after marginalization 
over both W and B. The dots indicate the regions containing 
90, 95, and 99 per cent of the probability, respectively 
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